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1 Introduction 



Recently LSND group [l[] (See also [fj]) reported that they have found candi- 
date events for z7 M — ► V e oscillation. Combing their result in Q with the data 
by the E776 group on the same channel — > u e and the reactor data by 
Bugey 0, the mass squared difference Am 2 seems to be in the range 

0.27eV 2 <Am 2 SND <2.3eV 2 . (1) 

On the other hand, it has been argued that the solar neutrino problem ||- 
U is solved if a set of the oscillation parameters of the two-flavor neutrino 



mixing satisfies one of the followings H]-[14]: 



(Am 2 , sin 2 26 Q ) 

(0(l(r lo eV 2 ), Oil)), (vacuum oscillation solution) 

(O(10- 5 eV 2 ), C(10" 2 )), (small - angle MSW solution) (2) 

(O(10- 5 eV 2 ), 0(1)) (large - angle MSW solution). 



Furthermore, it has been reported by the Kamiokande group | 15flfl6 l that 



their atmospheric neutrino data suggests neutrino oscillation with a set of 
parameters (Am 2 , sin 2 26) ~ (1.8 x 10~ 2 eV 2 , 1.0) for (v e <-> 2/ M ), (1.6 x 



10- 2 eV 2 , 1.0) for «-> u T ) (See also 0-0). 

The results of experiments on these anomalies have been given in the 
framework of two-flavor mixings in the original literatures, and have been 



analyzed from the viewpoint of the three flavor mixing by many people |2T 
||. It has been shown that strong constraints are obtained if two of the 
three anomalies are taken for granted, and it seems difficult |27| (See also 
pg|j ) to account for all these three anomalies within the three flavor mixing. 

Here we investigate the possibility in which one sterile neutrino as well as 
three known flavors of neutrinos are responsible for these anomalies. Some 



features of this possibility have been discussed in the past [p9|-|0^|. In this 
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paper we analyze in detail the mass squared differences and the mixing angles 
in this scenario. It turns out that if the number N v of light neutrinos is less 
than 4, then because of strong constraints by the solar and atmospheric 
neutrino observations, accelerator and reactor experiments including LSND 
as well as big bang nucleosynthesis, all the 6 mixing angles are strongly 
constrained. In this case 4x4 mixing matrix is effectively split into 2x2 
Ve *-> v s and <-> v T matrices, and only the small-angle MSW solution is 
allowed. If we assume that these neutrinos are of Majorana type, then the 
upper bound on (m„J from neutrinoless double j3 decay experiments suggests 
that these neutrinos are not heavy enough to explain all the hot dark matter 
components. If N v > 4, then all the solutions to the solar neutrino problem 
are allowed, hot dark matter can be accounted for by neutrinos and there may 
be a chance to observe neutrinoless double f3 decays in the future experiments. 

It has been pointed out that physics of supernova gives constraints on 
sterile neutrinos [[!5|][[36|, but these constraints apply only to Am 2 which is 
larger than the mass scale suggested by the LSND data, so we will not discuss 
this point in this paper. 

In section 2 we present our formalism on oscillations among four species 
of neutrinos. In section 3 we discuss constraints from reactor and accelerator 
experiments. In section 4 constraints by big bang nucleosynthesis are con- 
sidered. In section 5 we give constraints from neutrinoless double (3 decay 
experiments. In section 6 we examine the possibility that neutrinos are hot 
dark matter. In section 7 we discuss consequences in case of N v > 4. In 
section 8 we give our conclusions. 
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2 Oscillations Among Four Flavors 



Let us consider a model with three flavors of neutrinos v e , u^, u T and one 
sterile neutrino u s , which is singlet with respect to all the gauge groups in 
the standard model. We assume that (u e , u T , u s ) are related to the mass 
eigenstates (z/ 1; u 2 , z/ 3 , z/ 4 ) by the following unitary matrix: 



U 



( Vex U e2 U e3 U e4 \ 



C/ M i C/,,2 u,* u, 



1x2 V fj3 



'/jA 



U T \ U T 2 U T 3 U T 4 

V u sl u s2 u s3 u s4 j 

ia' A 3 V3i~f' Ag „VGi5' A i5 i 



-V3i/3A 8e -ioA 3 



(3) 



Here 



,7T 



i? 34 (- - ^34)i? 2 4(^4)e i5lA3 i? 23 (^ - ^ 23 )e 



181X3 



(4) 



is a 4 x 4 Kobayashi-Maskawa matrix for the lepton sector, 
sin%, 



COS @ijiSij 



Rjk{0) = exp (iT. 



jfct 



(5) 



is a 4 x 4 orthogonal matrix with 



3 k hm 



% (5 je 5, 



km 



(6) 



2A 3 = diag(l, -1,0,0), 2VS\ 8 = diag(l, 1,2,0), 2V6\ 15 = diag(l, 1, 1, -3) 
are diagonal elements of the sw(4) generators. 

The four phases a', j3', 7', 5' in front of Vkm in © can be absorbed by 
redefining the wave functions of charged leptons. If there is a Majorana mass 
term, however, the three phases a, ft, 7 cannot be absorbed |]37|]-[[4T|. Since 
these factors a, /3, 7 are cancelled in the probability of neutrino oscillations, 
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they do not affect the results of neutrino oscillations, but as we will see later, 
they do affect the effective mass of v e and the constraint on the masses of 
neutrinos from neutrinoless double f3 decay. 

We can assume without loss of generality that 

2 , 2 , 2 , 2 (n\ 

m 1 < m 2 < m 3 < m 4 . (7) 

Three mass scales Am| ~ e>(l(T 5 eV 2 ) or e>(l(T 10 eV 2 ), Am^ ~ 0(l(T 2 eV 2 ), 
^ m LSND ~ C(leV 2 ) seem to be necessary to explain the suppression of the 
7 Be solar neutrinos || , the zenith angle dependence of the Kamiokande multi- 
GeV data of atmospheric neutrinos [|16j, and the LSND data |J,@ so we will 
look for solutions of neutrino oscillations among four species where the mass 
squared differences are the three mass scales mentioned above. 

In the present case there are six possibilities, and they are classified into 
two categories. One is a case in which two mass eigenstates have different 
degenerate mass scales: 

(ia) 

2 „ 2 2 

i 2 < m 3 ~ m 4 



with (Am 2 !, Am 2 3 ) = (Am|, Am 2 tm 



Ami = Am 2 SND for j = 3, 4, k = 1, 2, (8) 



(ib) 



i 2 ~ ml <^ m\ ~ ml 



with (Am 2 1; Am 2 3 ) = (Am 2 tm , Am|) 

Am% = Am 2 SND for j = 3,4,k= 1, 2. (9) 



4 It can be shown that the probability of neutrino oscillation in vacuum in the framework 
of more than three flavors is reduced essentially to the two flavor case, once one assumes 
the mass hierarchy p3[)|p4||, so that the analysis of jlj indicates that the mass squared 
difference in our case should satisfy (0). 
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Another is a case in which three mass eigenstates have degenerate masses 
while one eigenstate has a mass far from others: 



na 



(iib) 



nc 



(iid) 



ml ~ ml ~ ml <C ml 



with (Am 2 21 ,Aml 2 ) = (Ami, Am atm) 

Ami = Am LSND for j = 1, 2, 3, (10) 



ml ~ wij ~ mj <C 777.4 



with (Am 2 21 ,Aml 2 ) = (Am^ tm , Am| 

— A om 2 f^v „' _ 1 9 

'4? 



Am* =Amj> SND for j = l,2,3, (11) 



777^ <C 7772 — m 3 — m 4 



with (Aml 2 , Am%) = (Am|, Am^J 

Am' = Am^ SND for j = 2, 3, 4, (12) 



777^ <C 7772 — m 3 — m 4 



with (Am^ 2 , Am^) = (Ara^ tm , Am|) 

Am^ = A777^ SND for j = 2, 3, 4. (13) 

The latter four possibilities (ii-a) - (ii-d) are excluded by reactor and 
accelerator experiments, as we will show below. The case (ib) can be treated 
in exactly the same manner as (ia) by changing the labels l<->3, 2^4 of the 
mass eigenstates, so in sections 3 and 4 we will discuss only the case (ia) for 
simplicity. 
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3 Reactor and Accelerator experiments 



Let us discuss constraints from reactor and accelerator experiments. The 
probability P{y a — > 2/3) of the transition v a — > up in vacuum is given by 



P(v n 



V 2 

2 ]T Im (UoiUZjUfcUf,,) sin (A^L) 



i<jr' 



l-4^|t/ Q 



It/, 



sin 



i<j 



2 / fc-EijL 
V 2 



(14) 
(15) 



where a = e, /i,r, s are the flavor indices, Ai?y = Ei — Ej = \jy 2 + mf — 
mf)/2E = AmjJ2E is the difference of the energy of the 



p 2 + m\ 



two mass eigenstates. (|i~4"|) and ([15]) are exact expressions in vacuum, and 
they are simplified if we assume the patterns (|8|) - ([Tj^) with mass hierarchy. 

3.1 Mass Pattern (i) 



Applying the formula ( |T4] ) and (|T5|) to the present case (^) with mass hier- 
archy, we obtain the following expressions: 

2 (AE 3l L^ 



P{v„ 



P(Vn 



3=3 



sin 



V 2 

„il +\U t *\ 2 )(U l 



(16) 



,3, 2 + |f44r)sin 2 (^)CIT) 



where the CP violating terms have been dropped out in fllBf ). From (|17|) the 
negative result of a disappearance experiment measured at distances Li and 
L 2 gives a condition: 

r > 4|f/ Ql | 2 |f/ a2 |Va(Am 2 , 1 ) 

|2 i Irr |2\/rr |2 , Irr |2\,j? /a™2 



+ 4(£/ Ql | 2 + |f/ Q2 | 2 )(f/ Q3 r + |^a4r)e/«(Am^) 
+ 4|f/ Q3 | 2 |^4| 2 e/ Q (Am 2 3 ). 



(18) 
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f a (Am 2 ) in ( |i8|) is denned by 



f a (Am 2 ) = - 



I ~(Am 2 L x \\ I 2 /Am 2 L 2 

( sin — — — ) — ( sin 1 



\ AE \ \ AE 

l/sin 2 2# Q (Am 2 ) for f a {Am 2 ) > 1 
const. (Am 2 ) 2 for (sin 2 (^j^)) < 1, 

(19) 



where 



sin 2 



/ Am 2 V 



\ AE 



r°° [i™** , , Srl ,^.dag(E,q) , / Am 2 L 
-n T dE dq e(q)FJE) p \ ,qj sin 2 ' 
) jo Jo dq 



N a p(L) 'Jo Jo rfg \ 4£ /' 

(20) 

N aP {L) = n T dE dqe(q)F a (E) p \ > q) , (21) 
jo jo ag 

F a (E) is the flux of neutrino z/ Q with energy i£, tit is the number of tar- 
get nucleons, e(q) is the detection efficiency function for charged leptons lp 
of energy q, d<7p(E,q)/dq is the differential cross section of the interaction 
v@X — > ipX', and e stands for the largest fraction of the appearance events 
allowed by a given confidence level. 

Let us first discuss the reactor experiment In this case we have 
\Am 2 21 L/AE\ < 1, so that (|18|) becomes 

~ > (|t/e3| 2 +|^e4| 2 )(l-|t/e3| 2 -|^4| 2 )/Bu g e y (Am 2 1 ) 

+ |f/ e 3r|f/ e 4r/Bu g c y (Am 2 3 ), (22) 

where /Bugey (Am^-J stands for f e (Am^) in (^) in case of the Bugey exper- 
iment M. From the (Am 2 , sin 2 29) plot in [|j] we have 

6 x 1(T 2 < / Bu gc y (Am^) < 25 for 3 x l(T 3 eV 2 < Am 2 3 < 9 x l(T 2 eV 2 
11 < / BuB (Am 2 ) <25 for 0.4eV 2 < Am 2 x < 2.5eV 2 . (23) 
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Since each term in (JIT) is semipositive definite, from ( p2[ ) and fl23|) we get 

1 



|£4 3 | 2 + l^4| 2 < 



l 



1 - 1//Bugey (Ami!) 



or 



\U e3 \ 2 + \U e4 \ 2 > 



l + Jl- 1//Bugey (Am 



31. 



(24) 



(25) 



On the other hand, we have the following relation between the probabilities 
P iye v e\ A{xj) and P^ {y e — > u e ; A(x)) for the solar neutrinos transi- 
tions in the two and four flavor mixings, respectively: 



p(2) 



v e -(l-\U e2l \ 2 -\U e4 \ 2 )A(x\ 



1 - \U K . 12 



'e3\ 



\u, 12 



+ \u e3 \ 2 + \u ( 



el 



(26) 



where A(x) stands for the electron density, and we have averaged over rapid 
oscillations sin 2 (Am 2 ; L/4i?) for all ^ (2, 1). Eq. fl26|) can be derived in 
a way similar to the case of the relation between P^ and P (3) [|2j(43|. To 
account for the suppression of the 7 Be neutrino among all the solar neutrinos, 
P^ 4 ) {y e — > v e \ A(x)) cannot be larger than | [pi], so the second possibility 
fl2"5| ) is excluded. Using the relations U e3 = 
can be derived from the expression @), we have 

1 

40 



'13' °14 



Si3Ci4e°" 2 , f/ e4 = Si 4 e 3 , which 



B«(A m y<- for 0.24eV 2 <Am 2 1 <2.5eV 2 (27) 



or 



#1 3) ^14 < 9°. 



(28) 



Next let us consider the disappearance experiment of — > f44|. In 
this case we have \ Am\ l L / AE\, |Am 2 3 L/4.E| < 1, so that 
1 



u, 



mi i 



CDHSW 



(Am 2 x ) 



(29) 
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where /cdhsw (Am^) stands for (Am^) in (|19|) for the CDHSW data 
P|. The allowed value for /cdhsw (Am^) is given by 



1 < /cdhsw (AmjJ < 50 for 0.24eV 2 < Ara^ < 2.5eV 2 
and therefore we obtain 



IM 2 



\u, 



< 



1 - V 1 ~ 1 //cdhsw (AmlO 



or 



> - 

2 



1 ~ 1 //cdhsw (Am^) 



(30) 



(31) 



(32) 



Furthermore, we have the constraints from the atmospheric neutrino data. 
If we consider the probability P{y^ — > v^) for the multi-GeV data by the 
Kamiokande group, then we get 



P(v» -> ^) 



l-4|^ 3 | 2 |^| 2 sm 2 

-2(|^ 3 | 2 + |^4| 2 )(1-|^3| 2 -|^4| 2 ) 5 



(33) 



where we have used \Am\ l L/4:E\ C 1. To have the zenith angle dependence 
of the multi-GeV data, it is necessary for P{p i1 — > Py) to deviate from unity 
to some extent. Using the constraints sf 3 , sf 4 < 1/40 in (|28| ) and the relations 



U, 



til 



"Cl2C23C24Sl3e %&2 + ^8x282^^ ~ C\2C\^S\/i l S2A& %&Z (34) 



U, 



-C23C24Sl2Sl 3 e ^ ~ C 12 C 2 4S 23^^ ~ C 13 S 12 S 14 S 2^ *** , 



(35) 



we have 



\U, 



\u, 



fj3\ 



,2 2 
'23 C 24 



< 



yfi - V/ 



CDHSW 



31/ 



(36) 



10 



Because of the condition (pE|), the probability P(u e — > v e ) is close to unity for 
the mass squared difference Am 2 ~ Am 2 tm , so the only possible source for the 
atmospheric neutrino anomaly in our scenario is deviation of the probability 
P{ v n ~^ v ij) from unity. Using the technique in [25], we find that the region 
of the parameters |t^ M 3| 2 and 2 allowed for the multi-GeV atmospheric 
neutrinos data of Kamiokande at 90% confidence level is 

0<# 23 <50 o , 25° < 6 24 < 55°, (37) 

where we have used the relations 

= c 13 c 23 c 24 - s 13 s u s 2i e 1 ^- 5 ^ (38) 

= Ci 4 S24- (39) 

From ([37]) we have s 23 c 2i < 0.48, and hence the possibility (|31~1) is excluded. 
Let us now consider the LSND data . In this case we have 

\Am 2 21 L/4:E\, \Aml 3 L/4:E\ < 1 (40) 

so we obtain 

e = Ae\U e3 U; 3 + U e4 U;,\ 2 f LSND (Am^) , (41) 



with 



/i— >e 

2 



f = l/sin 2 2# LSND (Am 2 ) for Am 2 > O.leV 
j ~ const. (Am 2 ) 2 for Am 2 < O.leV 2 , 

where sin 2 2^lsnd(A?tt. 2 ) is the value in the allowed region in the (Am 2 , sin 2 29) 
plot in |[], and 25 < /LSND(Am 2 1 ) < 580. Using the conditions on s 13 , s u , s 23 , 
s 2 4, we obtain 

1 2 

(si3 + s u ) /lsnd (AmjJ (43) 
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and hence 



9°>max(# 13 ,M>0.8 o , (44) 
where we have also shown the upper bound 

3.2 Mass Pattern (ii) 

We now show that the mass patterns (|T0| ) - ([Tj^) are excluded by reactor and 
accelerator experiments. The arguments against the four cases are exactly the 
same, so for simplicity we will consider the case flT0|), where Ami ~ A?til SND 
for j = 1, 2, 3. If the flight length L of neutrinos in the experiment satisfies 
AE41L ~ 0(1), then the oscillation probabilities are given by 

P(v a -+vp) * 4|t/ Q4 | 2 |^ 4 | 2 sin 2 (^|^) (45) 

12 f-, I tj \2\-2(^ E ^ L 



P(u a ^u a ) ~ l-4|[/ a4 r(l-|M 2 )sm 2 ^— f-j, (46) 

where CP violating terms have been dropped in P(v a vp). Exactly with 
the same arguments as in the previous subsection, the data by Bugey in this 
case implies 

1 



itu 2 < 2 

the CDHSW data leads to 

nu 2 < i 



1 - \A - 1 //Bugey (Am^) 



(47) 



1 - v^I-V/cdhsw (AmlJ , (48) 



and the LSND data to 

1 = 4|^ e4 | 2 |^ 4 | 2 /LSND (Am^) • (49) 
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So Am^ has to satisfy the following constraint: 



4|£/ e4 | 2 IM 2 



< 



/lsnd (Amlj) 



1-4/1 



/Bugey (Am^ 



/cdhsw (Am^i) 



(50) 



It turns out that the condition is not satisfied for any Am^ in the 
entire region 0.27eV 2 < Am^ < 2.3eV 2 . Hence the mass patterns (|iy) - (|1~3|) 
are excluded by reactor and accelerator experiments. 



4 Big Bang Nucleosynthesis 



Let us now discuss the constraints by cosmology It has been shown |46|] ]|7 
that the transitions v a — > u s (a = e, //, r) have to be suppressed strongly for 
big bang nucleosynthesis to be consistent with the standard scenario, if the 
number N u of light neutrinos is less than four.! Here we follow the argument 



of [4"5| ||46|| and give a rough estimate on the allowed region of the mixing 
angles without detailed numerical calculations. Except in section 7, we will 
assume N v < 4 in the following sections. We will show in section 7 that 
if N v > 4 we have much weaker constraints by big bang nucleosynthesis 
compared to the case for N v < 4. 

To get a rough idea on the magnitude of the effective mixing angle of 
neutrino oscillations, let us consider the temperature dependence of the dif- 
ference AE of the kinetic term and the potential term V |^5[ |46| [0 . 

Am 2 



AE 
V 



6.3T 

y / 2(7/90) 2 7r 6 (2 + cos 2 6 W )G F 

a?>w w 



rpis 



(51) 
(52) 



5 For recent discussions on N Vl see EqlfH and references therein. 
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From ([51]) and (|52| ) we have 



AE V13MeV/ \ eV\ 

In case of the two flavor neutrino oscillation with a mixing angle 9, the 
interaction rate of sterile neutrinos through neutrino oscillations is given by 
'FN F 



= sin 2 29 M sin 2 (££ coll ) ^, (54) 

V 12 J 2-flavor tl 

where 9m is the effective mixing angle in the presence of interactions defined 



by 



2 a — s ^ n2 ^ 

Sm 26 M = (cos29-V/AEy + sin 2 29' (55) 



H 



45 m 2 ! 



.1/2 



id 



(56) 



is the Hubble parameter, g*(T) = 10.75 is the effective degrees of freedom of 
particles for m e <T < m^, 

£ = - (A£ 2 + V 2 - 2VAE cos 29) 1/2 (57) 
is the eigenvalue of the mass matrix 

e i<T2e diag (E 1 , E 2 ) e" 1 ^ + diag ( V, 0) , (58) 



4oii = l/r„ = 1/C(T)G^T 5 (59) 

is the collision length of neutrino, T v is an interaction rate for active neutri- 
nos, C{T) ~ 0.5 for m e <T < m^, and 



H ~ (l.9Mev) ' ^ 
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Notice that there is no enhancement of oscillations due to the MSW mecha- 
PTCfl in the present case, because 8m < 9 follows from V/AE < in (|55|). 



msm 



Since we are interested in the mass squared differences Am^ ~ O(10 5 eV 2 ) 
or C(10- 10 eV 2 ), Am 2 3 ~ £>(l(T 2 eV 2 ), Am 2 31 ~ £>(leV 2 ), we have to con- 
sider situations where AE^ = Am 2 fc /6.3T becomes comparable to the abso- 
lute value |V| of the potential, and this condition is satisfied for the critical 
temperature T ~ 2MeV (or 0.3MeV), 8MeV, 15MeV, respectively. Off these 
regions of T, the problem becomes simpler. If T 3> 25MeV, \V/AE\ becomes 
very large and the effective mixing angle is small, so that we have 



— J < sin 29m t? 

tl J 2-flavor tl 



sin 2 26 



AE\ ( T 



V J U.9MeV, 
T \- 9 /Am 2 ^ 2 



for >25MeV. (61) 

On the other hand, if T -C lMeV, \V/AE\ becomes very small and this case 
is reduced to the oscillations in vacuum: 

^ < sin 2 2d ( T V<1 

H/2-fiavor ~ \1.9MeV7 

for < lMeV. (62) 

In the above extreme cases there is little transition u s — > v a (a = e,fi,r) 
and we have no problem with big bang nucleosynthesis. So in the follow- 
ing discussions we will discuss only the region of temperature m e <T <m^. 
Since AE£ co u/2 ~ 250\AE/V\ for any region of T, the term sin 2 (££ co u) in 
the formula of neutrino oscillations can be put to 1/2 after averaging over 
rapid oscillations, as long as the temperature under consideration satisfies 
\AE/V\ ~ 0(1). 
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To analyze oscillations among four species of neutrinos, it is necessary to 
diagonalize the mass matrix 

M = Udmg (Ex, E 2 , E 3 , E A ) U" 1 + diag (V, cV, cV, 0) , (63) 

where we have taken into account the fact that v e has both charged and neu- 
tral current interactions while and v T have only neutral ones for m e < t < m M , 
and 

c = cos 2 6 w /{2 + cos 2 9 W ) ~ 0.28 (64) 

is the ratio of neutral current interactions to the total contribution. Here 
we consider three cases where the potential term becomes comparable to the 
energy difference \V\ ~ AE jk = Am 2 jk /6.3T for (j, k) = (2, 1), (4, 3), (3, 1). 
This corresponds to the critical temperature T ~ 2MeV (or 0.3MeV), 8MeV, 
15MeV, respectively. From the constraint ([28|) by reactor experiments we 
know that U e 3 and U e 4 are small, so we take U e s = V eA = in what follows 
for simplicity. Since we assume N v < 4, sterile neutrinos should never have 
been in thermal equilibrium, and we demand that the interaction rate T Us 
of sterile neutrinos through neutrino oscillations be smaller than the Hubble 
parameter H: 

at=e,fi,T 

< H, (65) 



where 



VI + 4sin 2 ^ + 8sin 4 ^i v ; 



In (|6"5 ) both charged and neutral current interactions are taken into account 



for u e , while only neutral ones are included for v„ and v. 
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Let us first consider the situation where \V\ — AE21. For the vacuum 
oscillation solution, we can give the same argument and arrive at the same 
conclusion as the MSW ones. In what follows, therefore, we will discuss for 
simplicity the case of the MSW solutions, which imply a critical temperature 
T ~ 2MeV. Now if T ~ 2MeV, the mass matrix (^) is diagonalized as 

M ~ UU M [(Et + cV)U + diag (A 1; A 2 , AE ai , AE 41 )} U^U~\ (67) 

where 

Ai, A2 



1 



U M 
A 

with 



AE 21 =f - [AE 21 - 2AE 21 (1 - cA)V cos 26 Q 



1 2 



1/2 



\u. 



.s3 1 



It/, 



,s 1 1 



tan 20^ 



cA)^sin 2# G 







(68) 
(69) 
(70) 

(71) 



A£ 21 -(l-c4)Vcos2# ' 

In this case difference of any two eigenvalues of M. is large compared to the 
collision length 4oii, so the factor sin 2 (££ C oii) can be put to 1/2, and the 
probabilities are given by 

I— A 



where 
sin 2 2(^ + £ ) 



2 



sin z 2( 



; 12 



(72) 



A(l - A) {2 - i sin 2 2« + ^0)} + y sin2 2< ^> ( 73 ) 



sin 2 25 P 







1 - 2(1 - cA)Vcos29 Q /AE 21 + (1 - eA^/A^ ' 

(74) 
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and we have parametrized the matrix elements as 



U s3 

u si 



with 



(75) 



X = arg(C/ s3 D r s4 ), 
ijj = &ig(U s3 /U si ), 
4> = tan _1 |f/ S 4/f/ s3 |. 

Thus we have the interaction rate of sterile neutrinos 

T Ua = [p( Us ^ Ue ) + B {P{v a ^ + P{v 8 ^ v T )}]r v 



(76) 
(77) 
(78) 



B{2A(1 -A) + — sin 2 20} + {^ A B ^ "=~ 2 
for T ~ 2MeV. 



--^(l-A)}sin 2 2(^ + ; 



(79) 



If T ~ 8MeV (i.e., \V\ ~ AE 43 ), the mass matrix (|63|) is diagonalized as 
A< ~ UU M l(E 3 + cV)l 4 + diag (-AS31, A 3 , A 4 )] U^U- 1 , (80) 



where 



A3, A 4 
U M 



with 



A£ 43 ^ 



A£?„ + 2AE A3 cAV cos 2^ + c 2 A 2 V 2 



J 13 



2 

1 2 

e^ll 



tan 2^ = AF C ^t' 
34 A£ 34 + c4Vcos! 



1/2 



(81) 
(82) 



(83) 



All differences of two eigenvalues of Ai except AE 2 \ are large compared to 
the collision length £ coa , and we put all the factors sin 2 (££ C on) to 1/2 but set 
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sin 2 (Ai?2i-4oii) to zero. The probabilities are given by 



v.. 







A 1 



where 



sin 2 2(0 -Ci) 



2A{\ - A) + y sin 2 2( 



sin 2 20 



7 34 J' 



1 - 2cAV cos 20/A£ 43 + c 2 A 2 1/ 2 /A£| 3 ' 
The interaction rate of sterile neutrinos is given by 

A 2 



(84) 
(85) 

(86) 



r ~ 

1 Vs — 



B 



2A{1 - A) + ^-sin 2 2( 



rv 



(87) 



for T ~ 8MeV 

If T ~ 15MeV (i.e., \V\ ~ A.E31), the mass matrix (|63|) is diagonalized as 

M ~ UU M [ExU + diag (V, -AE 31 , A 5 , A 6 )] E^ET 1 , (88) 

where 



-AE31 + cv _ L 1 



17. 



,1/ 





2 


/ 













V 






A 



2 



-\ ( u sA 
{ -u s3 



{AE 2 3l - cV) 2 + AcAVAE 31 



1/2 



oi. 



rr* 
u sl 

JT* 

JT* 

u s3 

JT* 

u sA 



OL- 



^1 ^ \ 



TJ* 

JT* 

TJ* 



(89) 
• (90) 



In ([H|) a±,/3± satisfy the normalization conditions 

A|/3| 2 + (1-A)H 2 = 1, 
and x± = P±/a± are the two roots of a quadratic equation 
Ax 2 + (1-2A+ AE 31 /cV)x - (1 - A) = 0. 

19 



(91) 



(92) 



From @g) it turns out that (UU M )ei = 1, (UU M ) ej = for j = 2,3,4, 
(UUM)sj — for j = 3, 4, and the only combination which appears in the for- 
mula of probability with non-vanishing coefficients (UU m) S j(UU 'M)* k {UlJ m)%^ 
(UUM)ek is (j, k) = (3,4), and the difference of these two eigenvalues of M. 
is large compared to the collision length £ coU . Hence we again put all the 
factors sin 2 {£l co \\) to 1/2, and we obtain the probabilities 

P{y s ^v e ) ~ (93) 

P(y s -> Vp) + P{v a -> i/ T ) 

~ 2|a+| 2 |a_| 2 |l - A + x+A| 2 |l - A + a;_A| 2 (94) 
2,4(1 -A) 
1 + 2(1 - 2A)cV/AE 3 i + c 2 K 2 /A4 ' 
The interaction rate of sterile neutrinos is then given by 

r 2c2A ^ - A ) r 

1 + 2(1 - 2A)cV/AE 31 + c 2 V 2 /AEl v ' 

for T ~ 15MeV (96) 

Combining the results (|79|), fl8?D and (|96| ) for three cases we have the following 
formula 

2BA{1 - A) 



r / t 



if V1.8MeV 
+ : 



1 + 2(1 - 2A)cV/AE 3 i + c 2 \/ 2 /AE| 1 
5A 2 sin 2 20/2 



+ 



2cAV cos 20/A£ 43 + c 2 A 2 1/ 2 /A£f 3 
(1 - A)(l -5A) sin 2 2# /2 



1 - 2(1 - cA)V cos 26 & /AE 21 + (1 - cA^/AE^ 
which approximately holds for any temperature m e < T < mu i. As i 



(97) 



111 
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we look for the extremum of the expression T U J H with respect to T. The 



6 Notice that results (|79|), (|37|), (pq) are reproduced for each case in fl97|): the first term 
-► 2c 2 A(l - A), the second ->■ c 2 A 2 sin 2 20/2 for T ~ 2MeV where V/AE31, V/AE i3 -► 0, 
the first term -> 2c 2 A(l - A), the third - (T/2MeV)~ 12 -> for T - 8MeV 
where V/AE 31 -> 0, |F/AS 2 i| -> large, the second - (T/8MeV)- 12 -> 0, the third 
- (T/2MeV)- 12 -> for T - 15MeV where \V/AE 31 \, \V/AE ia \ -> large. 
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extrema of the three terms in (B7I) are attained at different values of T which 
are far apart, so we can discuss the extremum of the total quantity (|97|) 
by examining each term separately. Each term in ( |9"7| ) has the following 
maximum value 

15 x3 /A_2 \V2 



max (1st term) 



/ V eV 2 



3ff 



X 



A(l - A) [1/3 + ((1 - 2A)/3) 2 ] 1/2 - (1 - 2A)/3 



1/2 



2A) 2 /3 + (1 - 2A) [1/3 + ((1 - 2A)/3) 



2] 1/2 

(98) 



max (2nd term) 



7.7^ 3 

D 



Am 2 3 



1/2 



3A 3 / 2 B 



2 x 10- 2 eV 2 



-1/2 



sin 2 20 [1/3 + (cos 20/3) 2 ] 1/2 + cos 20/3 



1/2 



max (3rd term) 



1 + cos 2 20/3 + cos 20 [1 /3 + (cos 20/3) 

3(1- A)(l- BA) 



21 1/2 



(99) 



Am 2 2l 



1.9V 
L87 I 10" 5 eV 2 



1/2 



8(1 - cA) 1 / 2 



sin 2 26 P 



x 



[1/3 + (cos2# /3) 



21 1/2 



cos 29q /3 



1/2 



\2lV2 



.(100) 



1 - cos 2 26 G /3 + cos 26 Q [1/3 + (cos 26 Q /3)\ 
The term in the second line in ( |9~8|) is a monotonically increasing function in 
A, so the first term in (J9Tf) gives a necessary condition 

3 U/2 / A™2 \ -1/2 



.4 



< 



(— 



c 1/2 /Am 2 / 



<1.3 x 10" 



(101) 



3^3 V 15 J B \ eV 2 

where we have used the property that the term in the second line in fl98|) is 
approximated as y/3A/A as A — > and the lower bound Am^ >0.27eV 2 in 
([!]). It is straightforward to see that the term in the second line in fl9~9|) is 
less 4/3 for any value of cos 20, and it follows from ( pLOlj ) 



max (2nd term) 



< 



2 x 10" 



Am 2 ^ 



1/2 



2 x 10~ 2 eV 2 



<1, 



(102) 
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where we have used the combined Kamiokande results of sub-GeV [^5j and 
multi-GeV atmospheric neutrinos, which suggest 5xlO~ 3 eV 2 < Am^ < 3x 
10 _2 eV 2 . So the contribution coming from T ~ 8MeV never brings v s into 
thermal equilibrium, as long as max T (lst term) < 1 in (^) is satisfied. 
The term in the second line in (|100|) is approximately less than 0.54 for 
sin 2 2^0 < 0.9 which is satisfied by both the MSW solutions. Thus we have 

/ Am 2 \ 1/2 

max (3rd term) < 0.24 { j^p j > ( 103 ) 

so v s is not brought into thermal equilibrium in case of the MSW solutions, 
since Am^ <9 x 10 _5 eV 2 is satisfied by both the MSW solutions. Exactly 
by the same argument as above, we can show that the vacuum oscillation 
solution is allowed, as far as (|100|) is concerne dS. From (HD and fllOTD it 
follows 

|f/s 3 | 2 ,|^4| 2 <1.3 x 10- 2 (104) 
From (0) we have 0.33 < c 2 4 < 0.82, so that 

|[/ s4 | 2 ~c 2 4S 2 4 >0.33 4, (105) 
where we have used the relations 

U s3 = c 13 c u s 2 3e~ tSl - c 2 4Si3Sus u e t{S2 ~ 5s) - C13C23S24S34 (106) 

U s4 , = Ci 4 C 2 4S 3 4- (107) 

From ( |104j ) and (|105| ) we conclude 



s 2 M < 1.3 x IO-7O.33 ~ 3.8 x 10~ 2 , (108) 



7 The critical temperature for the vacuum oscillation solution is T ~ 0.3MeV, so the 
effective degrees g*(T) of freedom becomes smaller than 10.75. Nevertheless, it turns out 
that max T (3rd term)^<O(10 _1 ) x (Am| 1 /10 _5 eV 2 ) 1 / 2 is satisfied also for the vacuum 
solution and this maximum value is obviously much less than 1. 
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or 



From 



and 



#34<11 C 



, we have 



U s3 \ 2 ~ si, < 1.3 x 10" 



"28 



(109) 



(110) 



or 



9 23 <6°. 



Ill) 



To summarize, for N v < 4 the mixing matrix looks like 



KM 



( c s e e ^ 

£ £ C-atm ^atm 

c £ Satm Catm 

V -% c e e / 



(112) 



where we have defined 6 & = 12 , a tm = $24, e stands for a small number, and 
all the CP violating phases have been dropped out because |#i 3 |, |6*i 4 |, (#23 1, 1 #34 1 
<C 1. The basic reason that we have a strong constraint such as ( |112| ) is be- 
cause we have demanded that one of the mass scales is given by AmL SND , 
which implies a severe bound on T V JH from big bang nucleosynthesis. 

It was pointed out in [45] that v e <-> v s oscillation is potentially dangerous 
because it could change the density of v e at T <T U and therefore could in- 
crease the nucleon freeze-out temperature. To avoid this effect, it is necessary 
for the probability P{y e — > u s ) to satisfy 

/ in 75 \ 1/2 
P{v P -»■ u a ))) = ( , 7777 ^ ) > 0.93 (113) 



10.75 + 1(^-3) 



where we have substituted G 2 F — > G 2 F {1 + (1 — P{y e — > u s ))/2 in the equilib- 
rium condition 



G F T 5 ~ g*{T) 1/2 T 2 / 



(114) 
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in case of v e <-> v s oscillation, while we have used N v < 3.9 |48] and have put 
g*(T) = 10.75 -> 10.75 + (7/4) (JV„ - 3) < 10.75 + (7/4) x 0.9 on the right 
hand side of (|114p. As we have seen above, the only case where v e v s 



oscillation occurs is T ~ 2MeV (See ( \FA) ), and using (|113|) and the condition 
A <C 1 we have 

sin 2 2(6^ + ) < 0.26. (115) 

The numerical value here differs slightly from the one in [Q, but in either 
case we arrive at the conclusion that the large-angle MSW solution is ex- 
cluded, since ( |74"| ) exceeds 0.26 for some region of T for 0.6 < sin 2 2# <0.9. 
Note that this argument does not apply to the vacuum solution, because 
the collision length is much shorter than the oscillation length. We note in 
passing that the large-angle MSW solution for u e <-> v s is excluded without 
referring to big bang nucleosynthesis, by combining all the data of solar neu- 
trino experiments, the earth effect and Kamiokande day- night effect 1 12] HOI- 



It has been shown in the two flavor analysis Jl4| that the vacuum oscilla- 



tion solution for the channel v e «-> v s is excluded at 95 % confidence level if 
the standard solar model is taken for granted. The mixing matrix (|112|) tells 



us that our model in case of A^, < 4 is described approximately by two-flavor 
mixings v e <-> v s and u T , so we conclude that the vacuum oscillation 

solution is also excluded in our scenario. Hence we are left only with the 
small-angle MSW solution for N v < 4. This is consistent with a naive antic- 
ipation that the small- angle MSW solution is the most preferable scenario, 
because it fits the solar neutrino data best |TT|— fH|. 



5 Double f3 Decay 

Assuming that neutrinos are of Majorana type, let us now consider the im- 
plication to the neutrinoless double (3 decay experiments. To discuss neu- 
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trinoless double (3 decays, we have to consider the magnitudes of neutrino 
masses instead of mass squared differences. From our assumption on the 
mass hierarchy Am^, AmJ 3 <C Am| 1? Am^, Am^, Am^, we have 



mi ~ m 2 , m,3 ~ 772,4 



m\ + Am^. 



(116) 



From the form of the mixing matrix (112) the two pairs (u ej u s ) and (i^, v T ) 
have basically no mixing with each other, and we have 



(m Ue ,m Us ) ~mi, im v ^m VT ) ~ ^mf + Am| 

for mass pattern (ia). 

In case of the mass pattern @ , the mixing matrix looks like 



;ii7) 



( e e c Q s Q \ 

Catm ^atm C C 
"Satm C a tm ^ ^ 



;ii8) 



so that we have 



(m Ue , m Us ) ~ \Jm\ + Am| 1; {m v ^m Vn ) ~ mi 

for mass pattern (ib). 



(119) 



It is well known |J7|-[[IT| that the CP violating phases a, /3, 7 in (|3|) cannot 
be absorbed by redefinition of the neutrino fields. So the effective electron 
neutrino mass, which is the e-e component of the mass matrix, is given by 



V KM D diag (nij) D V^ M 



mi 



e ia c% + e~ ia s% 



mi 



for mass pattern (ia) 



m\ + Amli 



e ia cl + e 



-ia Q 2 



ml + Amli 



for mass pattern (ib), 



(120) 
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where D = e lv ^T Al se 4/3 As e ia T A 3 is a diagonal matrix and we have used 
( PT7D , ( P^l ) and the condition sin 2 2^ ~ 6 x 1(T 3 for the small-angle MSW 
solution. The present upper bound on (m Ue ) from neutrinoless double (3 
decay experiments is |5(| H51J1 



(m^) < 0.68eV. (121) 
From (p0|) , flmp , ([TTTD and ([mj) we obtain 



mi < 0.68eV for mass pattern (ia), (122) 



m x < V(0.68eV) - Am|i 

< 0.44eV for mass pattern (ib), (123) 

where in case of ( |123| ) we have used the lower bound in (|1|). Note that 
^ m LSND — ^ m ii < 0.46eV 2 has to be satisfied in ( |123| ) for mi to be real. 

6 Hot Dark Matter 

It has been argued that a model with cold+hot dark matter with 5eV < J2 a ^ Q 
< 7eV seems to be in agreement with observations, such as the anisotropy 
of the cosmic background radiations, correlation of galactic clusters, etc. 



[31][52|, and efforts have been made to introduce sterile neutrinos to account 



for hot dark matter as well as other anomalies p9|-||34|| . Here we examine 
the possibility that neutrinos could be hot dark matter while satisfying all 
the constraints that we have obtained in the previous sections. 

Since we assume that sterile neutrinos have never been in thermal equi- 
librium, only three components of the mass matrix ( |1 1 7| ) or ( |119| ) contribute 
to the mass density. Using (|122j) and (|123j), therefore, we have the following 
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total mass bound: 



mi + 2 Jm\ + Amfi < 4.0eV 
for mass pattern (ia) 



2mt + sjm\ + Amfi < 1.6eV 
for mass pattern (ib). 



(124) 



From ( |124[ ) we conclude that in either case neutrinos are not heavy enough 
to account for all the hot dark matter components. 

7 N v >4 



The possibility of N v > 4 has been proposed recently |kJ, so we will also 
discuss this case for the sake of completeness. If N u > 4, then sterile neu- 
trinos should have been in thermal equilibrium, and we cannot put a strong 
constraint like the one for N v < 4. The only condition which has to be satis- 
fied for sterile neutrinos to be in thermal equilibrium is A > 1.3 x 10~ 2 . We 
still have constraints from the reactor experiments which imply that #13, #14 
are small, but we can no longer say anything about the magnitude of #23, #34. 
The large-angle MSW solution is allowed in this long as v s is in 

thermal equilibrium. As we noted earlier, however, the large-angle MSW 
and vacuum oscillation solutions are excluded for v e <-> v s in the two fla- 
vor analysis |y| [ 13j |14 1 , so if U is very close to a direct sum of two mixing 



matrices with v e <-> u s and <-> v T channels then it contradicts with so- 
lar neutrino observations. The analysis of solar neutrino problem with four 
species of neutrinos would be extremely complicated!, and it is yet to be seen 



8 The explicit calculation for solar neutrino problem with three flavors has been per- 
formed recently in pa], assuming mass hierarchy. 
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under what conditions of mixing angles the large-angle MSW and vacuum 
oscillation solutions are allowed. 

Here, for simplicity, let us consider an extreme case in which the oscilla- 
tions take place mainly in the channels v e <-> v T and v s . This solution 
is consistent with all the solar neutrino observations |L2 1 [ 13| |14| , and u s is in 
thermal equilibrium because the large mixing angle of i/p v s suggested by 
the atmospheric neutrino anomaly JL6| gives T V JH which is larger than one 
j4j| . The masses of neutrinos in this case are given by 

(m Ue ,m UT ) ~ mi, (m v ^m Vs ) ~ ym 2 + Am^ 



for mass pattern (ia) 



{m Ue ,m Ur ) ~ ymf+Am^", (m Vft ,in Va ) ~ mi 
for mass pattern (ib). 



(125) 



(126) 



If the solar neutrino problem is solved by the small-angle MSW solution, 
then from (|122|) and ( |123| ) we have 



2mi + 2 J ml + Am 2 
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< 



4.7eV for mass pattern (ia) 
2.2eV for mass pattern (ib). 



(127) 



So in this case neutrinos cannot account for all the hot dark matter compo- 
nents. 

If we take the large-angle MSW solution, on the other hand, instead of 
fll20D we have 



miyl — sin a sin 29 Q > 0.32m! 



for mass pattern (ia) 



ml + Am^ a/1 - sin 2 a sin 2 26 Q > 0.32 Jm\ + Am 2 
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for mass pattern (ib), 



28 



(128) 



where we have used the condition 0.6 < sin 2 26 & < 0.9 for the large-angle 
MSW solution. From (TM and (TM we have 



mi < 2.2eV 



2mi + 2y m 2 + Am^ < 9.6eV for mass pattern (ia), 

(129) 

m 1 < 2.1eV 



2m x + 2ym 2 + Am^ < 8.5eV for mass pattern (ib), 

(130) 

so that neutrinos could explain all the hot dark matter. Conversely, if all the 
hot dark matter components are neutrinos in this case, then it follows 



2m! + 2y/m? + Am 2 31 > 5eV, (131) 
mi > 0.8eV, (132) 
0.2eV for mass pattern (ia) 



0.4eV for mass pattern (ib), 

(133) 



so we will be able to observe neutrinoless double (3 decays in the future 
experiments. 

If the solar neutrino problem is solved by the vacuum oscillation solution, 
we have 



mi\/l — sin a sin 29 Q > 0, 

for mass pattern (ia) 



mf + Am 2 ! y 1 — sin 2 a sin 2 29 Q > 
for mass pattern (ib), 

(134) 
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where we have used the condition 0.6 < sin 2 29 P) < 1 for the vacuum oscil- 
lation solution. If a and 26 Q are both very close to tt/2, then mi can be 
arbitrarily large without contradicting with the bound ( |121| ) from the neu- 
trinoless double (3 decay experiments. In this case, therefore, neutrinos could 
account for all the hot dark matter components, but there is no guaran- 
tee that we will be able to observe neutrinoless double (3 decays in future 
experiments. 

8 Conclusions 

In this paper we have performed a detailed analysis of the constraints on a 
model with three active and one sterile neutrinos, using the data of reactor 
and accelerator experiments, the solar and the atmospheric neutrino obser- 
vations, and big bang nucleosynthesis. The mass pattern where three masses 
are degenerated is found to be inconsistent with reactor and accelerator ex- 
periments. If N v < 4, then all the mixing angles are severely constrained 
and the mixing matrix is effectively split into 2x2 matrices with channels 
v e *-* v s an d <-> v T . In this case the large- angle MSW and vacuum oscil- 
lation solutions are excluded. Because of the constraints from neutrinoless 
double f3 decay experiments, neutrinos cannot explain all the hot dark mat- 
ter components. For N v > 4, we get fewer conditions on the mixing angles, 
and all the solutions to the solar neutrino problem are allowed. If we take ei- 
ther the large-angle MSW solution or the vacuum oscillation one, then these 
neutrinos can account for all the hot dark matter components. In case of 
the large-angle MSW solution for N u > 4, we will be able to observe neutri- 
noless double (3 decays in near future. It is hoped that combined results of 
super-Kamiokande and SNO experiments can tell us about the existence of 
sterile neutrinos, the type of the solution to the solar neutrino problem, and 
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the mixing of sterile neutrinos p3| . 
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